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Abstract 

An improved treatment of neutrino oscillations follows when neutrino mass eigenstates and their 
associated recoiling particle states are entangled. When the neutrino and its recoil partner are 
treated as a single entangled quantum state the conservation of energy and momentum in neutrino 
oscillations can be assured, even in a plane- wave treatment. The oscillation wavelength between the 
neutrino and its associated recoil partner emerges as the fundamental periodicity in this analysis. 
However, the experimentally determined oscillation wavelength for neutrinos detected at some 
distance from a known source region is still in all cases given by the standard expression in current 
use. 
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INTRODUCTION 



Neutrinos produced in a weak charged- current decay are prepared in a state of speci- 
fied flavor because they are accompanied by a charged lepton with that flavor. As is now 
well established, such flavor states are in fact a linear superposition of mass eigenstates of 
mixed flavor. The mass differences lead to a phase slip between the quantum waves of each 
mass component, which can be experimentally detected as a change in flavor at some place 
spatially separated from the point of decay. 

Neutrino oscillations, a phenomenon recognized in Pontecorvo's farsighted prediction in 
1957 [lj, can be treated in a simple and straightforward way leading to a correct expression 
for the wavelength in terms of the neutrino energy and the difference between the squares of 
the masses of the eigenstates. But a propagating neutrino that consists of 2 or more mass 
eigenstates cannot simultaneously have a well defined energy and a well defined momentum. 
How Lorentz transformations should be made then also depends on which is chosen. These 
difficulties are well known - a thorough summary has been given recently by Akhmedov 
and Smirnov [2]. Generally the solution is to introduce a wave-packet approach, which is 
necessary for other reasons in any case, and treat both the energy and the momentum as 
somewhat undefined. Another strategy is to use quantum field theory, treating the neutrino 
as a propagator and thereby avoiding detailed consideration of an internal line, but also 
rendering it difficult to obtain oscillations at all [31 0]. 

Such difficulties do not arise when the neutrino and its recoil are treated in a unified 
way as entangled-state components. The role of entanglement in neutrino and neutral kaon 
oscillations has been considered by others [31 I5HT2]. who, with some exceptions [TOj IT2] . 
show that the standard formulae are obtained. In our analysis, an oscillation wavelength 
between the neutrino and its associated recoiling partner emerges. When a source decay 
region is then localized, that oscillation wavelength can be related to a laboratory oscillation 
wavelength for the neutrino distance from the source, which can be shown [T3| ITS] to be the 
same as the standard wavelength. 

The basic principles are best seen in two-body decays, such as pion decay, or electron 
capture. Two-body decays are interesting for another reason: - the kinematics that produce 
neutrino flavor oscillations produce, under restrictive conditions (of which one is that the 
neutrino be observed), an 'oscillation' of the recoil particle or nucleus. 
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CALCULATION 



Consider a two-body decay at rest, for example 7 Be + e — > 7 Li + v e . The final state is 
a two-particle state comprising the neutrino and the recoiling 7 Li ion: 

\Ru e ) ~ \R) \p e )p- 1 exp i{pp-Et). (1) 

It has the form of a product wave function of the internal degrees of freedom of the neutrino 
and recoil particle and a function of their relative motion, a spherical wave. The energy E 
is the energy of the final state, which is the same as the energy of the initial state. The 
momentum of each particle is p and their separation is p. This is the basic form of the two- 
body wave function that will be recast in Lorentz invariant form below. The neutrino and 
recoil particle are entangled in a single quantum state. A key feature of this entanglement 
is the appearance of a single conserved energy E and a single time coordinate t. 

The (electron) neutrino is a linear combination of mass eigenstates with masses m«, 



:> = E^k>, (2) 



where the U e i are elements of the Maki-Nakagawa-Sakata-Pontecorvo mixing matrix [15J. 
In the rest frame of 7 Be (mass M), decaying into an electron neutrino u e and a recoiling 
particle (henceforth simply 'the recoil', in this case 7 Li) with mass m, the neutrino (recoil) 
energies are [15] 

M 2 - (+)m 2 + (-)m 2 
E im = — . (3) 

The sum of neutrino and recoil energy is always M in the rest frame (the masses include 

the appropriate numbers of atomic electrons). The momenta are equal and opposite for the 

recoil and neutrino eigenstate, 

([M 2 - (m + m,) 2 ][M 2 - (m - m,) 2 ]) 1 / 2 
IPil - 2M ' O 

A complete description of the final state must include the wave functions of both the neu- 
trino, which may have a general flavor projection and the recoil. The wave function can 
be written schematically in terms of the four-momentum P and the spacetime coordinates 
X of the recoil, and the corresponding quantities p and x for the neutrino: 

\Rv) ^Y.UIMP.X))^^)),. (5) 



Kinematics requires the subscript % on the recoil wavefunction. The recoil and the neutrino 
move as free particles, asymptotically in plane-wave states. 



w r (P)e- itrP ^w s (p)e- 



(6) 



The w are spinors, e = +1(— 1) identifies positive-energy (negative-energy) states, and P,p 
are the 3-momenta of the recoil and neutrino respectively [16]. The particles move along 
the x-axis. The spin degrees of freedom and negative-energy states are irrelevant for the 
present, kinematic, considerations. In the rest frame of the parent, P = —p. 



i 

= \R)e-^Y.u:M) 



iE R t—ip x X ^—iEvt+ipxX 



Jp x (-X+x) 



(7) 



where E a = M is the energy of the parent. The final state is seen to have the desired form 
of a product wave function of the internal degrees of freedom of the neutrino and recoil, and 
a function of their relative motion, as in Eq. [TJ 

Viewed from a frame moving at velocity — /3 along the x-axis, the momenta have the 
transformed values P' x = •yf3En — jp x and p' x = •y(3E u + 7p x , the space coordinates become 
X' = 7/3t + 7X and x' = j{3t + jx, and the time coordinates become t' v — jt + r yf3x and 
t' R = yt + ^PX. The energies become E' R = ^E R — •yf3p x and E' v = r yE u + •y(3p x . The primed 
quantities represent the quantities in the laboratory frame when the parent is boosted to 
a velocity (5. Since the arguments of the exponentials in Eq. [6] are Lorentz invariants, one 
may equivalently write 

W rs/ = Ei U* ei \R) h) x 

w r^ e -ie r P'^' w s^ e -ie 3 p'^ 



\Rvy^EiU* ei \R)Wi) 



e -iE R t R +iP x X e -iE v t v +ipxX 



If the experimental situation calls for observing the space and/or time coordinates of the 
particles by detecting them, then the laboratory spacetime coordinates are of interest. The 
phase 

\-E R t R + P X X - E v t v + p x x]' 
= -(E R + E v )t + T l p x {-X' + x'). 



The energy En + E v = E a is the energy of the final state in the rest frame, which is 
the same as the energy of the initial state, and therefore carries no index i. The internal 
coordinates of the recoil are independent of % as well. When the neutrino is detected, a 
measurement is made of its relative probability density at a specific spacetime location, 
assuring that the spacetime coordinates of the point of detection are independent of i. The 
state detected is: 

\Ru)' = e~ lEat \R) J2 Ki \n) e^~ lp ^- x ' +x,) . (8) 

i 

The leading phase factor that depends on the energy is unobservable. As a result, the 
interference effects of neutrino oscillations arise solely from the different momenta in the 
components in the final state. (In this we concur with |17j . but see also [H]). Since p Xj i 7^ p x j, 
phase differences exist between the components at the point of detection. The electron flavor 
admixture for two admixed components i,j oscillates with a wavelength A^- given by 

([M 2 - (m + mi ) 2 ][M 2 - (m - m,;) 2 ]) 1 / 2 2vr 
([M 2 - (m + m,) 2 ][M 2 - (m - m^ 2 }) 1 ' 2 

Dropping terms of order m\ and writing m 2 — m 2 = Am^ , 

X'ij Airhc 7M M 2 — m 2 ' 1 ; 



where the missing units have been made explicit. All quantities in Eq. 10 except 7 are 
Lorentz invariants. In terms of the 'neutrino beam energy' E' Q , which may without approx- 
imation be defined as the energy that massless neutrinos would have in the laboratory, 

1 _ Am 2 . 1 + p M 2 + m 2 

X'ij ~ 4nhc E' 2M 2 ' { ' 



Equations 10 and 11 define a wavelength in laboratory coordinates for the separation — X'+x' 
between the neutrino and the recoil [TH] . 

In the standard expression [T5] for neutrino oscillation, the survival probability oscillates 



with a wavelength A£ ■ given by 



1 Am 2 - 



It appears the two results are different, but Eq. 11 defines an oscillation length between 



neutrino and recoil, whereas Eq. [12] is the oscillation length relative to the decay point. 



Experimentally, the latter situation is the usual case and so it will be desirable to recast 



Eq. 11 into an equation for the wavelength of oscillations in x' alone. This can be done by 
taking into account the fact that X' and x' are not independent variables, but are related 
by kinematics once an origin has been defined. We may conveniently choose the origin in 
spacetime as the point at which the decay occurs. In effect, this defines a phase relationship 
between the plane waves such that the flavor is electron at the origin. In doing so, we 
appear to run afoul of the uncertainty principle because plane waves have precise momenta 
and indefinite positions. Replacing the single plane-wave states i,j with distributions of 
plane- wave states having closely similar momenta (by imposing a small boost distribution on 
the parent), one creates correlated wave-packets. These superpositions of plane- wave states 
nevertheless retain the general properties of interference identified for the single plane-wave 
states, provided the width of the momentum distribution is not so large that the interference 
pattern is made indistinct at the point of detection. 

The oscillation probability may be found from Eq.[8j noting that \R) describes the internal 
coordinates of R and thus is independent of i: 



P{Ve^Vp) = \(Rup\Ru e )\ 
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(13) 



which may be expanded in the usual way |15j . In the oscillation probability, the i—j interfer- 
ence term depends on (—X f + x') through the phase factor exp [«7 _1 (pa:,t — Px,j)(— X' + x')]. 
It is noteworthy that in the entangled-state picture, oscillations arise even in the plane-wave 
basis because states containing different neutrino eigenmasses nevertheless have the same 
total momentum (zero in the rest frame) and the same total energy. 

To find the oscillation probability in the more customary coordinate x' , when the decay 
occurred at or near the origin, the variable X' must be expressed through kinematics in 
terms of x'. In the rest frame of the parent, the 3-momenta are equal and opposite, and the 
distances are proportional to the speeds: 

2M 2 

-X+x~x— (14) 

M 2 + m 2 K ' 

where the neutrino masses have been neglected (their contribution turns out to be higher 



order than mf [II] [19]). Referring to Eq. 11 , one sees that for (3 = 0, this replacement yields 



a wavelength in x that is identical to the standard one, as Vogel has shown [T3] . 



This equivalence extends, in fact, to all frames, as shown by Kayser [14J. A general 
argument will be presented elsewhere [19J, but there is an interesting special case for which 
it is clear. If the parent is boosted by a velocity that leaves the recoil at rest in the laboratory, 
then the wavelength for neutrino oscillations relative to the decay point is given directly by 



Eq. 11 The speed of the recoil in the rest frame (again neglecting the neutrino masses) \/3r\ 



is given by 



2M 2 

1 + = XrTT^' (15) 

M z + m z 



and it is evident that for this particular boost, Eq. 11 once again gives a result equivalent to 



the standard result, Eq. 12 That the standard result arises from an entangled-state ansatz, 
without disentanglement, has been shown also by [7HH]. 

The entanglement with the recoil enforces a specific relationship between the momenta 
of the oscillating neutrino components, and does not permit an arbitrary choice (usually 
either 'beam energy' or 'beam momentum') about which to expand in order to derive the 
oscillation phase. The total energy, rather than neutrino eigenmass energy, contributes to 
the state phase, removing it from a physical role. An illustration of the relationship is shown 
in Fig. [Tj 

The group velocity of the mixed state in this formalism is equal to the classical particle 
velocity: 

AE _ |p| 
A|p| ~~ ~E' 

The importance of this condition for a physically reasonable description has been emphasized 
(e.g. Levy [20]); here it is derived. 

The form of Eq. [8] shows that oscillation effects are manifested in the separation between 
the neutrino and the recoil. Whether 'recoil oscillations' are observed depends on whether 
and how the associated neutrino is detected. In a realistic situation where the location of 
the decay is known to some accuracy, the recoil will always be found at a distance from a 
detected neutrino given by kinematics and it will show no oscillations whether the neutrino 
is detected in coincidence with the recoil or not. In a less practical case where the decay 
point is unknown, the recoil-neutrino distance will show periodic oscillations in space (in 
the rest frame of the parent). Our conclusions match those of Dolgov et al. (6]. Consider 
a toy model of 2-neutrino mixing with equal amplitudes of electron and mu flavor in each 



m. 




FIG. 1. Illustration of the relationship between energy and momentum for two neutrino components 
of mass mi an d mi- The two lines with slope +1 are the Einstein equation for particles with 
those masses. The standard formula for neutrino oscillations may be obtained by selecting a 
common average momentum (line a) or energy (line b) and expanding to derive an oscillation 
phase corresponding to the energy or momentum difference, in time or distance coordinates. Both 
happen to give the same result, because the slope is +1 (to second order in . In the entangled-state 
picture, the momentum difference between two neutrino mass eigenstates is given by the difference 
in 3-momenta at the intersections with line c. The expression displayed shows the modification of 
Am 2 when the phase in the neutrino-recoil system is calculated (in a frame other than the rest 



frame of the parent, an additional factor 1 + /3 enters - see Eq. 11). Despite this constraint, as 
Cohen et al. [5j show, the phase for neutrino oscillations relative to the origin is always given by 
the standard expression, for any momentum/energy combinations allowed by the eigenmasses. 

eigenmass, i.e., 

\Ru) ~ \R) [h) + e^teM-fe.iX-^'+z') i^)] ( 16 ) 

omitting leading phase factors and constants. The neutrino may be detected as an electron 
neutrino, for which the phase is 0, 2ir, etc.; as a muon neutrino, for which the phase is 7r, 37r, 
etc.; by a neutral-current process, for which the phase is irrelevant; or not detected at all, 
for which the phase is again irrelevant. Detecting the neutrino by a charged- current process 
enforces a certain phase relationship, and if the recoil is detected in coincidence with the 
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neutrino, its probability will show spatial oscillations as well because the separation — X'+x' 
between recoil and neutrino may be increased by an integer number of wavelengths to recover 
the same state. On the other hand, when the neutrino is detected via the neutral current, 
or is not detected, then the phase is irrelevant, detection of the recoil is unaffected by it, 
and the probability of detecting the recoil does not depend on X'. A situation with some 
similarities was noted by Smirnov and Zatsepin [21] for neutrinos emitted in Z° decay, where 
both particles are neutrinos. 

CONCLUSIONS 

We turn to a summary of our conclusions. Taken in isolation, a propagating neutrino 
state in a superposition of mass eigenstates may be assigned a well-defined energy or a well- 
defined momentum, but not both. Equally troubling, the appropriate way to make Lorentz 
transformations depends on the choice. Ignoring those difficulties, it is possible to derive 
in a simple way the standard form for neutrino oscillations (see [6] for a summary of such 
derivations). On the other hand, when the neutrino and its associated recoil are treated as 
a single state whose components are entangled the above objections drop away. Oscillations 
emerge as entirely a kinematic effect observable as a spatial pattern in the wavefunction of 
the two-body entangled state, and Lorentz transformations are unambiguous. 

In a previous version of this paper, we found the expression for the wavelength of os- 
cillations in the parent-recoil separation coordinate and erroneously applied that to the 
oscillation probability of the neutrino detected at some distance from a source whose loca- 
tion was approximately known. The wavelength in the neutrino-recoil-separation coordinate 
differs from the standard form of the wavelength in the neutrino-origin coordinate for neu- 
trino oscillations. In fact, while the two wavelengths are different, they are related in a 
way that makes them completely consistent [13, 14]. Entanglement in this case therefore 
produces no readily observable physical consequences. In principle, a measurement of the 
recoil momentum would permit an unambiguous prediction of the eigenmass detected in a 
neutrino detector, but such precision has not yet been realized, nor are neutrino coincidence 
experiments simple to do. It is often objected that entanglement does not apply to situ- 
ations where a 'measurement' is made 'before' the neutrino is detected, because the recoil 
hits something. We would argue that the central requirement in a valid quantum-mechanical 
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analysis is to deal with an isolated system, and it is that which forces entanglement on us 
in the first place. It may be formally necessary to enlarge the system to include subsequent 
interactions of the recoil. The conservation of energy that removes the energy from the phase 
still applies. The momentum relationship between neutrino components will apply only un- 
der the condition that the surroundings play no role until the neutrino and recoil are beyond 
the range of their interaction. There are clearly cases where that is not a good assumption, 
for example, the known influence of a solid crystalline medium on beta decay [22J . Moreover, 
a complex system that encompasses the recoil and has internal degrees of freedom would 
technically invalidate the assumption that \R) has no dependence on the index i, although 
the range argument could be invoked to rescue it. It is therefore fortunate, perhaps, that 
entanglement makes no alteration to the standard form of neutrino oscillations. 

The momenta and energies of components of a propagating neutrino state have a defined 
relationship that is fixed by the properties of the associated recoil particle. This relationship 
leads to a group velocity that corresponds to the velocity of the neutrino. The interference 
of the components emerges even in a plane-wave description and is due to the differences 
in momenta. The oscillation probability has a spatially periodic dependence on the sepa- 
ration of the neutrino and its associated recoil. The wavelength depends on the masses of 
the parent and recoil, and on the Lorentz boost of the parent. However, transformed into 
coordinates for the position of the neutrino relative to the decay point, the oscillation prob- 
ability corresponds to the standard form and depends on neither the masses of the parent 
and recoil nor on the Lorentz boost. 
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